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Abstract 

The QCD evolution equations in Modihed Leading Log Approxima- 
tion for the factorial moments of the multiplicity distribution in quark 
and gluon jets are numerically solved with initial conditions at threshold 
by fully taking into account the energy conservation law. After applying 
Local Parton Hadron Duality as hadronization prescription, a consistent 
quantitative description of available experimental data for factorial cu- 
mulants and factorial moments of arbitrary order and for their ratio 
both in quark and gluon jets and in e~^e~ annihilation is achieved. 
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1 Introduction 



The analytical perturbative approach to multiparticle production in jets 
is widely used to describe inclusive observables in e+e~ annihilation. In this 
framework, a parton cascade is evolved down to small scales of a few hun- 
dred MeV for the transverse momentum cutoff Qo, and partonic predictions 
are directly compared to hadronic observales, according to the notion of lo- 
cal parton hadron duality (LPHD) 0. The description of the parton cas- 
cade is given within the so-called modified leading logarithmic approximation 
(MLLA), which takes into account the leading double logarithmic and next 
to leading single logarithmic terms; in this approximation, QCD coherence is 
included by means of angular ordering, and the running of the coupling as at 
one-loop level and energy conservation effects are included as well. 

This approach has been applied to a large variety of experimental data and 
found to be rather successful (see for a recent review). In particular, by 
considering predictions which fully take into account the effects of the energy 
conservation law, a quantitative description of available experimental data for 
single particle observables, like for instance the momentum spectrum |^ and 
both the jet and the particle multiplicity 0, has been achieved. However, 
in the case of observables related to genuine multiparticle correlations, like 
factorial cumulants and factorial moments of the multiplicity distributions, 
the available theoretical predictions |^ , which include only partially the energy 
conservation effects, do not quantitatively reproduce the experimental data^]. 

In this letter we extend the results of for the average multiplicity and we 
present a numerical solution of the complete MLLA evolution equation of QCD 
for factorial moments of any order with the full inclusion of energy conservation 
effects. In this way, we can quantitatively describe within the purely pertur- 
bative approach recent experimental data on multiparticle correlations both 
in single quark and gluon jets and in e~^e~ annihilation, thus lending further 
support to the analytical perturbative approach to multiparticle production. 

2 The theoretical framework 

In order to study genuine multiparticle correlations in a quark or gluon jet, it 
is convenient to study the properties of the (unnormalized) factorial moments 
of order q, given in terms of the multiplicity distribution P„ by 
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and the (unnormalized) factorial cumulants of order g, related to the factorial 
moments by a cluster expansion: 
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Factorial moments and factorial cumulants are integrals of g-particle inclusive 
distributions and g-particle correlations functions, respectively. Normalized 
observables are obtained by using the average multiplicity n = F^^^ = K^^^: 

Also the ratio of factorial cumulants over factorial moments 

H, = KJF, (4) 

has been widely studied. 

One usually also considers the generating functions of the multiparton final 
states in quark and gluon jets, Za{u), from which the factorial moments can 
be obtained via a Taylor expansion around u = 1 (here a = q,g for quark and 
gluon jet respectively): 

and similarly for the cumulant moments 

The evolution equations for the generating functions Za{u) have been de- 
rived within a probabilistic description of the parton splitting processes A — > 
B + C, with the inclusion of angular ordering, energy conservation and the 
running coupling at the one-loop order Ijl],^,^ (we explicitly show the energy 
dependence and drop the label u everywhere for clarity): 

r dZ,{r]) _ /-i-c asih). 
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dz'^^[^liy{z){Z,{r^ + lnz)Z,(r/ + ln(l - ;.)) - Z,{r^)]+ 
nf^9,{z){Z,{r^ + ln^)^,(r7 + ln(l - z)) - Zg{r^)]] (7) 

= dz'^^%,{z){Z,{^ + In z)Z,{^ + ln(l - z)) - Z^} 
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Here the evolution variable r/ = In ^ is related to the jet virtuality k = Q sin ® , 
where Q = 2E is the hard scale involved in the process, and E is in this case 
the jet energy, G denotes the maximum angle between the outgoing partons 
B and C and Qq is the infrared cutoff in transverse momentum. 



The splitting functions ^ab for parton splittings A ^ B [|10[ are taken 
with normalization as in here = {l — z)^gg{z) can replace ^^gg{z) 

thanks to symmetry properties of the kernel |^ . Nc and n f denote the number 

of colours and flavours. The coupling is given by = 27r/(61n(^x/A)) 

with b = (llA^'c — 2?7./)/3. as{k±)~^ evolves with k± with a smooth treatment 
of the heavy quark thresholds at the effective mass m* = ~ l.lSmj 

as in ip. The boundaries of the integrals are determined by the lower cutoff 
in the transverse momentum measure k± defined according to the Durham jet 



finder algorithm |]Tl| kj_ = mm(z, 1 — z)k > Q^. 
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By taking the derivatives now in eq. (0) according to eq. (^, one then 
obtains a system of coupled evolution equations for the unnormalized factorial 
moments for a single quark or gluon jet: 

' dXgiri) _ f^-^c as{k 



dz'^^[^liy{z){Ug{r^ + \nz) +Ug{v + ln(l - z)) - Ar,(r/)}+ 
nf^g^{z){U^{r] + \nz)+ U^{r] + ln(l - 2;)) - Ug{r])}] 



dr] 



dr] 



f dz%^%g{z){Afg{r] + \nz)+ Af,{r] + ln(l - z)) - Afg{r])} 



dr] Jzc 27r 

{ E fl VfH^ + \nz)F^^-"'\r] + ln(l - z)) - (r^)}+ 

q 



m=0 



m=0 



dr] Jzc 2n 



m=0 



19 

m 



{E ]^'''^\r] + \^z)Fj^'^-^\r] + \ri{l- z)) - F^'^\7])} 



3 



Since Zc = e ^ < ^, the system of differential equations is defined for 
?7 > In 2 only. Its initial conditions for < r/ < In 2 read 



p{q) = = forg > 1 (10) 

No analytical solution of the complete equations (||) with the boundary condi- 
tions ([T0|) has been so far obtained. A numerical solution of the two coupled 
equations for the average multiplicity has been given in For high or- 
der factorial moments, both solutions within high energy approximations at 



leading fT^ and next-to-leading order |jT3[ and a solution of yet higher order 



with partially takes into account the energy conservation law^ have been ob- 
tained. Notice that these solutions do not fulfill the absolute normalization at 
threshold, where they reach unphysical negative values for the moments. 

Here we numerically solve the system of equations (P) with boundary con- 
ditions at threshold (^^) for the first 20 factorial moments. 



2.1 From one to two hemispheres 

We have so far considered particle production inside a single quark or gluon 
jet. In this framework, a single quark jet is defined in an inclusive sense, i.e., 
it experimentally corresponds to a single hemisphere in a e~^e~ annihilation 
event. At low cms energies this factorization of the two hemispheres does 
not hold anymore and nonlogarithmic corrections become important. At large 
cms energies, however, the behaviour of a whole e'^e~ annihilation event is 
controlled by the generating function 

Ztyjo-hemiu) = [Zgiu)]^ ^ In Ztyjo-herniu) = '2\n[Zg{u)] (11) 

It is then easy to relate the moments for the whole event to the moments for 
one hemisphere that one can calculate according to the evolution equation (H). 
One obtains indeed for the factorial cumulants: 

'^two-hem = , ^two-hem ~ '^K'^^ (12) 

and for normalized cumulants 

two— hem 2'"^ V / 

Factorial moments can then be calculated from the factorial cumulants by 
inverting eq. (|]). For instance, for the second order factorial moment, one 
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Quark jet 


Gluon jet 




Expi 


Theory 


Exp0 


Theory 


2 


1.0820 ±0.0006 ±0.0046 


1.080 


1.023 ±0.008 ±0.011 


1.026 


3 


1.275 ±0.002 ±0.017 


1.265 


1.071 ±0.026 ±0.034 


1.078 


4 


1.627 ±0.005 ±0.042 


1.600 


1.146 ±0.059 ±0.074 


1.157 


5 


2.274 ±0.014 ±0.093 


2.168 


1.25 ±0.11 ±0.13 


1.268 



Table 1: Experimental data for the normalized factorial moments of order 
g = 2, . . . , 5 in single quark and gluon jets of 45.6 and 41.8 GeV respectively P 
are compared with our theoretical predictions with parameters given in 
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The simple treatment of the whole e^e~ event in terms of a superposition of two 
independent single quark jets have been improved for the average multiplicity 
by explicitly including the 0{as) matrix element for e~^e~ —>■ 3 partons[0,0. 
In this way, an improved description of data in the low cms energy region 
could be achieved 0]. Since we are here mainly interested in the high energy 
regime of the moments, where the correction is negligible, we will use in the 
following the aforementioned simple ansatz. 



3 Phenomenology of high order correlations 

3.1 Correlations in single quark and gluon jets 

Both the evolution equation (|^ and eq. (|]) refer to parton production inside 
a single quark or gluon jet in a fully inclusive configuration. This is a very 
important point for the phenomenological application of these predictions. It 
is indeed not possible to directly compare our predictions with data extracted 
in different configurations, like for instance data on single quark and gluon jets 
from 3-jet events with the Mercedes configuration, where topology has been 



shown to play an important role [16,17| . However, it is possible to study a 



particular subsample of 3-jet events, with a very energetic gluon-jet practically 
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0.07431 


0.01922 


0.00135 


-0.00211 
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0.000330 


-0.000197 


-0.000125 


-0.000024 
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-0.00018 


0.00056 


0.00052 


0.00003 


-0.00041 


-0.00043 


piq) 
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0.000033 


0.000004 


-0.000005 


-0.000005 


0.000007 


-0.000009 
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14 


15 


16 


17 


18 


piq) 
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0.00004 


0.00060 


0.00066 


-0.00023 


-0.00155 


-0.00169 


piq) 

9 


-0.000015 


-0.000018 


0.000012 


0.000001 


-0.000024 


0.000067 



Table 2: Predictions for the ratios Hg in single quark and gluon jets of 45.6 
and 41.8 GeV respectively with parameters given in (|TBp. 



back-to-back to the quark-antiquark pair, where the experimental configura- 
tion is comparable to the theoretical one|]TB|. This analysis has been performed 



by the OPAL Collaboration both for the average multiplicity||T9[| and for higher 
order factorial moments|Q for quark jets of 45.6 GeV and gluon jets of 41.8 
GeV. Table (|l|) compares the experimental data for the first four factorial 
moments inside a single quark and gluon jet with the theoretical predictions 
obtained by numerical solution of eq. with the boundary conditions (]T0|). 
The free parameters of the theory, i.e., the infrared cutoff Qo, the QCD scale 
A and the normalization parameter Kaii, have been kept equal to the values 
fixed in to describe the jet and particle multiplicities: 

K^^i = 1, A = 500 MeV Qo = 0.507 MeV (15) 

A very good quantitative agreement between data and the predictions of the 
complete theory is achieved, contrary to previous comparisons with approxi- 
mate predictions 0. 

To be more explicit. Fig. 1 compares the data for factorial cumulants 
of order g = 2 up to 5 in single quark jets of 45.6 GeV and single gluon 
jets of 41.8 GeV with the predictions computed in this paper and with the 
predictions of . It can be seen that the full inclusion of energy conservation 
is very important for these observables and that, by including this effect, the 
perturbative approach can quantitatively describe experimental data in this 
configuration. 
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It is also possible to obtain predictions for high order moments in single 
quark and gluon jets. Predictions for the ratio Hq inside single quark and 
gluon jets up to order 18 are shown in Table 2. It is particularly interesting to 
notice that oscillations of the Hq moments are up to two orders of magnitude 
larger in quark jets than in gluon jets. 



3.2 Correlations in e+e annihilation 

Fig. 2 shows the dependence on F = hl{^/s/QQ) of the second order nor- 
malized factorial moment F^'^\ Experimental datapOl are compared with the 



leading and next-to- leading predictions and with our predictions, where 
the result for the full event has been obtained from the prediction for a single 
quark-jet according to eq. (|T3|) in the following way: 

= 1 + ^ (16) 

Also in this case, our improved description is in quantitative agreement with 
the experimental data and reaches an accuracy comparable to the widely used 
Monte Carlo generators. 

In order to investigate whether our predictions give a good description of 
correlations of yet higher order as well, we have also compared them with 
the data for the ratio of factorial cumulants over factorial moments, Hq, as a 
function of the order q. It has been shown in |2l|] that an important part of the 



oscillatory behaviour of the Hq moments actually comes from the superposition 
of two different samples of two-jet and multi-jets events, an effect which is not 
taken into account in our simple treatment of correlations between the two 
hemispheres according to ([T3|) . At the moment, we cannot then aim at a 
description of data for the whole e^e~ annihilation events [^, which requires 
an improved treatment of the composition of the two hemispheres, but we 
should look at data for two-jet events only. We have therefore compared in 
Fig. 3 our predictions with data for the Hq moments extracted as in ||21| from 
experimental data[^ on multiplicity distributions in three samples of two-jet 
events, defined according to the JADE jet finder algorithm with three different 
values of the jet resolution parameter ycut- One can easily see that the theory 
can describe the gross features of the data for correlations of arbitrary order 
as well. In this case, a more quantitative comparison should not be pursued, 
since data have been extracted by using a different jet finder algorithm than 
the one used in the theory. 
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4 Conclusions 



We have numerically solved the QCD evolution equations in MLLA for factorial 
moments of any order in quark and gluon jets, by fully taking into account the 
effects of energy conservation law. By keeping the same parameters previously 
fixed in order to describe the mean jet and particle multiplicities, a quantitative 
description of existing data on factorial moments, factorial cumulants and their 
ratio both in single quark and gluon jets and in the two-jet sample of e^e~ 
annihilation has been achieved. The analytical perturbative approach, based 
on perturbative QCD description of parton evolution assuming Local Parton 
Hadron Duality, is then shown to be quantitatively successful in describing 
genuine multiparticle correlations inside jets. 
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Figure Captions 



Fig. 1: a Data on the first 4 normalized factorial cumulants for single quark- 
jets of 45.6 GeV as measured by the OPAL collaboration (diamonds) com- 
pared with predictions by Dremin et al. p (triangle) and with our predictions 
(square); b: same as in a, but for single gluon jets of 41.8 GeV. 

Fig. 2: Second order normalized factorial moment F^'^^ as a function oi Y = 
\n{^/s/2QQ, with Qq = 0.507 GeV; comparison of experimental data|^ with 
leading order predictions |12|] (dashed line), next-to-leading order predictions ||13|| 
(dotted line) and our prediction (solid line). 

Fig. 3: Ratio Hg of factorial cumulants over factorial moments as a function of 



the order q; data for two-jet events as extracted in from experimental data 



on MD's measured by DELPHI Collaboration with three different values 
of jet resolution parameters within the JADE algorithm (?/c = 0.01 (diamond), 
0.02 (triangle) and 0.04 (square)) are compared with our predictions (solid 
line) . 
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Figure 1: a 
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Figure 1: b 
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Figure 2: 



13 



0.0005 p 
^ 0.0004 - 
0.0003 - 
0.0002 - 
0.0001 - 
- 

-0.0001 - 
-0.0002 - 
-0.0003 - 
-0.0004 - 
-0.0005 - 




_L 



8 



10 



<> 

-A— I 

T 

z\ 4^ 



12 14 16 18 



Figure 3: 
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